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Abstract 

This paper concerns the inverse problem of retrieving a stationary potential for 
the Schrodinger evolution equation in a bounded domain of with Dirichlet data 
and discontinuous principal coefficient a{x) from a single time-dependent Neumann 
boundary measurement. We consider that the discontinuity of a is located on a 
simple closed hypcr-surface called the interface, and a is constant in each one of the 
interior and exterior domains with respect to this interface. We prove uniqueness 
and lipschitz stability for this inverse problem under certain convexity hypothesis 
on the geometry of the interior domain and on the sign of the jump of a at the 
interface. The proof is based on a global Carlcman inequality for the Schrodinger 
equation with discontinuous coefficients, result also interesting by itself. 
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1 Introduction 

The method of Carleman estimates was introduced in the field of inverse problems by 
BuKHGEiM and Klibanov in reference [8] (see also [7] and [E]). The first known 
results concern uniqueness of inverse problems. Then, one of the first stability result for 
a multidimensional inverse problem, dealing with an hyperbolic equation, can be read 
in [25] and is based on a modification of an idea of [8]. 

Carleman estimates techniques are presented in [18| for standard coefficient inverse 
problems for both linear and nonlinear partial differential equations; one can also read 
in this book the construction of globally convergent numerical methods for coefficient 
inverse problems and some concrete applied problems in geophysics, medical imaging 
and computational time reversal. 
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It is possible to obtain local Lipschitz stability around the single known solution, 
provided that this solution is regular enough and contains enough information (see |17| 
and [IS])- Actually, many of the results using the same strategy we can refer to concern 
the wave equation. A complete list is too long to be given here but to cite some of 
them, related to the same kind of inverse problems of determining a potential and also 
using local or global Carleman estimates, see [23] and [27] for a Dirichlet boundary data 
and a Neumann measurement and [11] for a Neumann boundary data and a Dirichlet 
measurement. 

Recently, global Carleman estimates and applications to one-measurement inverse 
problems were obtained in the case of variable but still regular coefficients, see [13] for 
the isotropic case, and [H] and [3] for the anisotropic case. It is interesting to note 
that these authors require a bound on the gradient of the coefficients, so that the idea 
of approximating discontinuous coefficients by smooth ones is not useful. Nevertheless, 
uniqueness and Lipschitz stability are obtained in [1] for the inverse problem of retriev- 
ing a stationary potential for the wave equation with Dirichlet data and discontinuous 
principal coefficient from a single time-dependent Neumann boundary measurement. 

One can also note that a global Carleman estimate has been obtained [11] for the 
heat equation with discontinuous coefficients. That work was initially motivated by the 
study of the exact null controllability of the semilinear heat equation, but the estimate 
has been recently used (see [4] ) to prove local Lipschitz stability for a one measurement 
inverse problem. In this field, one should also read the recent works [5], [6] and [23]. 
Up to our knowledge, the result of determination of a time independent potential in 
Schrodinger evolution equation with discontinuous principal coefficient from a single 
time dependent measurement on the boundary is new. Concerning the simpler case of 
a "classical" Schrodinger equation (with a = 1), one can have a look at [2], where the 
Carleman estimate and the proof of the stability of the same inverse problem are maybe 
easier to read and the philosophy is the same. For the same equation, one can find in [22] 
a method with weight functions satisfying a relaxed pseudoconvexity condition, which 
allows to prove Carleman inequalities with less restrictive boundary observations than 
in [2]. The authors of ^0] deal with Carleman estimates for the Schodinger equation 
with variable (but regular) principal coefficient and applications to controllability. Let 
us notice that in the different context of Cauchy problem, V. ISAKOV in [15] uses local 
Carleman estimates for the Schrodinger equation to prove uniqueness of the solution. 
Finaly, for the Schrodinger operator idt + div(cV) in an unbounded strip in M^, reference 
[9] gives a stability result for the diffusion coefficient c in with only one observation 
in an unbounded domain. One will see in the proof of our main tool (an appropriate 
Carleman estimate) that it is based on the same strong pseudoconvexity condition (-^4) 
for the weight tp. 
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1.1 Statement of the problem and main results. 

Let T > and let C (A^ > 2) be a bounded domain with C^-boundary dn. 
Throughout this paper, we use the fohowing notations : 

^ / dv ^ A 

\dxi ' " ' ' dxN J ' ^ dxf 

, dv „ d'^v 

u G denotes the unit outward normal vector to 557, 

— = Vv.v is the normal derivative. 
ov 



We will work with the following Schrodinger equation : 

' iy'{x,t) + d\Y{a{x)Vy{x,t)) +p{x)y{x,t) = 0, x £ t £ (0, T) 
< y{x,t) = h{x,t), x£dn,te{0,T) (1) 

, y{x,0) = yo{x), X eQ. 

We consider in this paper the inverse problem of the determination of the coefficient 
p of the lower order term in Schrodinger equation ([1]) from a single time dependent 
observation of Neumann data ^ on the boundary. 

The major novelty of this paper is that we deal with a Schrodinger equation in a 
bounded domain of with discontinuous principal coefficient. Indeed, let Q and l^i 
be two open subsets of with smooth boundaries T and Fi. We choose Q,i simply 
connected and such that Jli C and we set = Thus, we have di}2 = T U Fi 

and we also set: 

fll X £ 0,1 
02 X G ^2 



a{x) 

with aj > for j = 1, 2. 



Considering equation ([T]), we know that for each p G L°°{^1), yo £ L^{Q) and h £ 
L^{r X (0,r)), there exists a unique weak solution y such that 

y G Ci[0,T];H~\n))nH-\0,T;L\n)). 

The proof is based on a transposition method, as one can read in [21] . Let us also notice 
that the regularity of y also implies — G H~'^{0,T; H~2 (T^). 

We will prove the well-posedness of the inverse problem consisting in retrieving the 
potential p involved in equation ([1]), knowing the flux (the normal derivative) of the 
solution y{p) of (pQ) on the boundary. It means that we will prove uniqueness and 
stability of the nonlinear inverse problem characterized by the nonlinear application 



I dy 
P\a ' — > 02 — 



(2) 

rx(o,T) 
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^ rx(o,T) 



in suitable norms ? 



Indeed, we will only give a local answer about the determination of p, working first 
on a linearized version of the problem, as shown is Section 3. Assuming that p G L°° is a 
given function, we are concerned with the stability around p. That is to say p and 
are known while q is unknown. We can also add that uniqueness is a direct consequence 
of stability but historically, uniqueness results were obtain first (see [8]) and stability 
was proved using for instance compactness- uniqueness arguments as in |27j . 

In this work we introduce a Carleman weight whose spatial part is similar to the one 
of the weight function constructed in [1] for the two-dimensional case. We prove a new 
Carleman estimate for the Shrodinger equation (see Theorem [9]) under the hypothesis of 
strong convexity -also called uniform convexity- for the interface (roughly speaking, it 
means that their curvatures are uniformly bounded from below by a positive constant; 
see Definition [6|), and some sign for the jump of the main coefficient. The following 
result, based on this Carleman estimate, states the stability of the inverse problem. 

Theorem 1. Assume that Qi is strongly convex and ai > 02 > 0. Let U be a bounded 
subset ofL°°{n), p e L~(f]) and r > 0. // 

2/0 € H^{Q.) is real valued (or pure imaginary) and if 



q GU, where y{p) is the solution of equation (1) with potential p. 

The main idea is that the nonlinear inverse problem is reduced to some perturbed 
inverse problem which will be solved with the help of a Carleman estimate. In order to 
obtain such an estimate, we first rewrite ([1]) as a system of two Schrodinger equations 
with constant coefficients and solutions yi and y2, coupled with transmission conditions. 
We then construct a Carleman inequality on each domain with nonzero boundary values 
on the interface. Next, we gather all the terms to construct a global Carleman inequality 
for the transmission problem. The main point is to look carefully at the interaction of 
yi and y2 on the common boundary Fi. 

Notice that we state hypothesis for a function which guarantee that it would be a 
suitable weight function for a Carleman estimate with the only requirement that the 



yoix)\ > r > 0, a.e. in Cl 




//i(0,T;L2(r)) 
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discontinuities of a are located on Ti. We shall only construct an explicit weight func- 
tion for the case of a discontinuous coefficient which is constant on each subdomain (i.e. 
oi and 02 constants). However, we could also construct a weight function for variable 
coefficients ai{x) and a2{x) such that their traces at the interface are constant, under 
additional assumptions of boundedness of Vaj similar to those appearing in [T3] (in 
order that the corresponding weight function would satisfy hypothesis (H^) and (-^4) in 
section [2]) . 

This article is organized as follows. Section 2 is devoted to the proof of an appropriate 
global Carleman inequality and Section 3 concerns the proof of the Lipschitz stability 
of our the inverse problem. 



2 A global Carleman estimate 

In this step, we will show a global Carleman estimate concerning a function v = v{x,t) 
equals to zero on dQ x {—T,T) and solution of a Schrodinger equation with a bounded 
potential q = q{x). We set the following notations : 

Q = nx{-T,T) 17o = f^iU02 

T = dn Ti = dni n 
s = rx(-r,r) Si = rix(-T,r) 

and if n is a function defined in Q, for Uj we will mean its restriction to the set flj, for 
each j = 1, 2. 

The main hypothesis for the Carleman estimate is the existence of a weight function 
ij; = ij;{x) defined on such that, on the one hand it is pseudo-convex with respect 
to the Schrodinger operator in each one of the two sub-domains and Q.2-, and on the 
other hand it has a convenient behavior at the interface Fi. 



Indeed, we will first suppose that G C^(0) verifies the natural transmission con- 
ditions: 

-01 = ^"2 on Ti 

ax- h a2T^ — = on V\. 

We will also suppose the following behavior at the interface 

ip{x) = cte for all x G Fi, (^1) 

dijji 811)2 , . 

In the interior 0,q we will need that 

>6>0 {H3) 

and that 3 e > such that 

2Z)2^(e, + 2a^X\Vij ■ ei' - aVa ■ > elCl' (^^4) 
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e C", where 



2 , I d ( dip 



■"J / / l<i j<iV 

Finally, it will be useful to consider weight functions satisfying (-ffs) and (-^4) except 
in a neighborhood of a point. In this case we will need two weight functions ip^ and if)'^^ 
each one satisfying {H^) and (^^4) in U 1^2 \ Bf,{xi) and f^i U ^^2 \ Bi,{x2) respectively, 
(with £ small enough) such that 

for each j, k G {1, 2} with j ^ k. 

Summarizing, we set the following 

Definition 2. 1. Let U C 0, be an open set such that Fi C U and let tp G C'^{U\Ti). 
We say that ip is a transmission weight function for equation ([I]) inU if it satisfies 
the conditions (TV), (Hi) and {H2) on the interface Vi, and hypothesis {H3) and 
{Ha) in U. 

2. Let 'p^ and V'^ be two functions in C^{Q.i U 0.2)- We say that {ip^ ,'ip'^) is an e-pair 
of transmission weight functions for ([I]) if there exist xi,X2 G and e > such 
that for each k = 1,2 the function ip^ is a transmission weight function for ([7]) in 
fio \ Bf,{xk) and the hypothesis {H5) is fulfilled. 

Given ip, for s > 0, A > we define on Q = Q x {—T,T) the following functions: 
6{x,t) = — -— and ip{x,t) 



{T-t){T + t) ' ' iT-t)iT + t) 



where a > \\e^^\\L°°{n)- 



We also define the space 



f dv ^ 

Z= lv£ L^{-T,T;H^{n)) : Lv £ L^{Q), — G L^{^) and v satisfies (Tr) \ , 

introduce the following norm in Z 

\\w\\s,\^ = s^X'^ f f e^\w\^dxdt + s\ f f e\Vw\^dxdt (3) 

and for || • ||s,a,i/),c/ we will mean the above terms defined in the set U C il.. 
We finally set 

Lv = iv' + div(a(x)Vu) + qv, 
V = e^^w 

and 

Pw = e-"^L(e"^w). 
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Hence we have 

Pw = iw' + isip'w + div(oVii;) + 2saVip.Vw 
+ sw div(aV(/9) + s'^a\\7(p\'^w + qw 
= PiW + P2W + qw 

where we denoted 

Piw = iw' + diY{aVw) + s'^a\V(p\'^w, 
P2W = isif'w + 2saS/ Lp.Vw + s div(aV(/3)u). 

Our main result is the following 

Theorem 3. Suppose there exists for some e > an e-pair of transmission weight 
functions belonging to C^{Qi U $72)- Let 9^, ip^ and w^ be the corresponding 

functions defined for tp^ as we did before. We also define 

= G r X (-r,r) : v^''{x,t) • u{x) > 0} 

Then there exists C > 0, sq > and Aq > such that 
<C (iP^'iw") _ +s\ 



+ 



w 



k=l 



dw^ 



d 



V 



for all V £ Z , X > Aq and s > sq. 



(4) 



2.1 Formal computations 

We have 

\Pw — qwl"^ dxdt 



I Iq 



Piw\^dxdt+ I I \P2wfdxdt 



+2Re J J PiwP2wdxdt, 
where z is the conjugate of z and Re{z) its real part. 

As 7; G L'^{-T,T;H^{n)) and v' e L'^{-T,T; H^^{n)) (because Lv G L'^{Q)), we have 
V G C{[-T,T];L^{a)) and w E C{[-T,T]; L'^{n)) with w{x,±T) = 0. 

We look for lower bounds for 



i?e yy P1WP2W dxdt = {Piw,P2w)j^2 
3 

We set {Piw,P2w)^2 = ^i,j-> where /jj is the integral of the product of the ith- 

term of Piw and the jth-term of P2W. The properties of w and some integrations by 
parts allow to write the following equalities. 
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To begin with, we have 

III = Re J J iw' {—isLp'w) dxdt = —- J J ip"\w\'^ dxdt. 

Applying the identity Im{z) = —Im{z) for z = 2s\ j j OaVip-Vww' dxdt we obtain: 

Ii2 = Re jj iw' {2saV(p ■ VW) dxdt 

= sX Im 9{ div(aV'0) + Xa\V'il^\'^)ww' dxdt 



sX Im / / aO'wVij) ■ Vw dxdt 



Q 

If— dijj 
+sX I I luaOw'—— dadt 
JJj] du 

We also have 

/i3 = Re j j iw' s div{a'V(p)w dxdt 



'21 



sXIm / / 6{ div(aVV') + Xa\Vij\^)ww dxdt, 



Re JJ d\Y{aVw){—isLp'w)dxdt 

sX Im 1 1 aO'wVip ■ Vw dxdt + s Im [ [ ip'wa^^ dadt, 



and 



JJQ J JT, 

I22 = R^ JJ dW{aVw){2saV(p • Vw) dxdt 

= - JJ Oa\Vw\^{ div{aVij) + Xa\V'ip\'^)dxdt 

- sX 9a\Vw\^Va-V'iljdxdt+ 2sX^ ii da^\Vip.Vw\'^ dxdt 



+ 2sX ReJJ eDl{ip){Vw,Vw)dxdt 



IQ 

2sX JJ 9a^{Vip-Vw)^dadt + sX JJ ea^\Vw\'^^ dadt. 
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Using integrations by parts we obtain 

-^23 = Re J J div(aVt(;)s div{aVip)w dxdt 

= sxjj \Vw\'^9{a div(aVV') + AjaV^p) dxdt 



^11 \w\'^ div{eaV div{aV'tp))dxdt 
Q 



2 



+ ^ 1 1 \w\'^OaV diviaVip) ■ udadt 



s 



2 



div(a^V(a|V^Ap)) + div(a^ div(aV'0)VV')) fixdt 



+— II \w\'^9a{dw{aV^p)VTp + V{a\VTp\'^)) ■ i^dadt 



''^^ f f I |2 1- / 2/1IV7 ; |2v7 ; N J Jj. i '^"^^ /" /" I 1 2 2/)i w / |2 '^V' 



— yy \w\' d\-v{a'e\V^'Vij)dxdt + — j j \wya'e\V^\'-^dadt 



-XsRe J J wOai div {aV ip) + XlVipl^ a)— da dt. 
and we obviously have 



I31 = Re 1 1 s'^a\Vip\'^w{—isip'w) dxdt = 0, 



I32 = Re J J s^a\V(f\^w{2saV(f ■Vw)dxdt 

= s^^^ jj \w\'^O^a{\Vij\'^ div{aVi;) + 2aD'^{il^){Vij,Vi;)) dxdt 
+ s^A^^y" \wf0^a\Vij\'^Va-Vijdxdt 

.3\4 ff |„.,|2/)3„2|v7„;,|4 . \3„3 ff I„„|2„2|v7„/,|2/)3'9V' 



+ 3s^A* // \w\^e^a^\Vil]\'dxdt-ys-' j j \w\'a'\Vxlj\^e-'-^dadt, 



and 



/33 = Re J J s a\V(p\ w{s dw{aV(f)w) dxdt 
= - s^X^ JJ \w\'^9^a\Vi;\'^ div{aVi;)dxdt 



Then we have Re J J P1WP2W dxdt = F{w) + GCVw) + J + Xi where we define 

F{w) = 2s^A^^y" \w\'^e^a'^\Vijfdxdt, (5) 

G{Vw) = 2s\^ [ [ ^a^l W • Vwp dxdt + 2sXRe [ [ 9DI^P{Vw, Vw) dxdt 
J Jq J Jq 

-s^jj \Vw\'^eaVa-V'ipdxdt, (6) 
J as the sum of all boundary integrals 



dip 

J = sXIm I I aOw'w-— dadt 



du 

sIm [ [ (p'wa—— dadt 
JJt, ou 

-2sXRe / / a^eVi; ■ Vw— dadt 
JJt. di^ 

+sX 1 1 a^e\Vw\'^^dadt 



s 



dv 



dvj 

-sXRe I I w div( a9'V'il>)a—- dadt 

„3\3 / / ^2/)3|„,,|2|v7„/,|2 ^V' 



-s-'X-' // a'9^\w\^\Vi;\^^dadt 
J Jt. ou 

ff 2/11 i2ivT ,|2^V' 



s 



[[ a^e\w\^\ViP\^^dadt 
J Jt. ou 



+— // a^e\w\'^y{\V^'^)-udadt 



s 



2 



H / / ae\w\'^ dw(aVil))—- dadt 

2 JJ^ du 

+yJJ a9\w\'^V{div{aVi{;)) ■ i^dadt 

and Xi as the sum of all the remaining integrals in 

Moreover, if [/ C is an open set, we will write Fu{w) to denote the sum of integrals 
from the definition of F{w) taken in the set U, and the same for G, Xi... 

Noticing that 

• 2sX Im j j 9'awVij.Vwdxdt 

<sX [[ a^{9')^Vij.Vw\'^dxdt+ sX [ [ {9')^w\^ dxdt, 
JJq j Jq 

• a G W^'°°{n) and ^ G C^(n) 

• 1^1 < C9'\ \9'\ < G9^ and|/| < G9^ on {-T,T) x Q, C = C{T) > 0. 
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it is then easy to prove, from simple calculations, that the "negligible" terms Xi indeed 
satisfy 



l^il <Cs\ o^^l W • Vw\^ + CsA^ / / d\w\^ + Cs^X^ / / ^VP- (7) 

J Jq JJQ J Jq 

2.2 Proof of the Carleman estimate 

In this part of the paper, we prove Theorem [31 We apply the above computations in 
each one of the domains Jli and 0,2 and we sum up all the terms. Since the interface Fi 
has null -measure, we get an estimate in all the set fi, plus the boundary terms from 
di}, and from the interface itself, where appear terms coming from both l^i and 
Given the hypothesis we have assumed, we prove in the following propositions that we 
can deal with all this terms. In the sequel, C denotes a generic constant, depending on 
T and Q. 

2.2.1 The interior 

Recall the norm || • \\s,x,^l) defined in ([3]) and F{w), GCVw) defined in 

Proposition 4. Suppose that U C il. is a open set and i{j satisfies (H^) and {H/^ in 
Oi U O2 \ U . Then there exist 7 > , C G M, sq and Aq such that for all v & Z , 

F{w) + G{Vw) + Xi> 7||t(;||s,A,t/- - C\\w\\sXi>,u 
\ls > So and VA > Aq. 
Proof : 

First, merely by the fact that ip E (7^(0), we have 

\Fu{w)\ + \Gu{yw)\ + \Xi^u\ < C\\w\\s,x,^^u (8) 

for all V £ Z. 

Now, from ([7]) and ip satisfying (^^3) we get that for s and A large enough, 

JJQ JJQ 

Hence, If satisfies (-ffs) and (-^4) in 17* = 17i U \ we get that VA > Aq, s > sq 
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and V G Z, 



> Fn.H + Gn.(Vu;)-|Xi,n. 



> FnAw) + Gn,(yw)-sX^ [ [ ^a^lW • V 

J-T Jn, 

|w;|2eV|VV'|^ 



T 

-T Jn* 



rT r rT 



+ 2sXRe [ [ 6Dli;{Vw,Vw) - sX [ [ \Vw\^eaVa ■ Vil^ 



> s'X' 



[ \w\^e^a^\Vij\^ + esX [ [ e\Vw\ 



> 7\\w\\sX^,n*- 
Prom ([8]) and Q we get the desired result. 



(9) 



2.2.2 The boundary 

By definition we have = on the exterior boundary S for each v £ Z. Therefore, 
dw 

Vit; E = -7^1^ and if we choose the legitimate notation J = + Je. , we get here 
ou 



dw 



2sXRe // a^eVilj ■Vw—dadt + sX II a^e\Vw\^^dadt 



-sX 



> -sX 



> -sX 



dw 
ou 



di/j 



dw 



du 



> -sXC 



dw ^ 



dv 
dadt 

dv 

dadt 

du 



dw 
du 



dv 



dadt 



dv 



dadt 



where we have denoted S+ = G F x [—T,T) : V'ip{x,t) ■ ^{x) > 0}. 



(10) 



2.2.3 The interface 
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We compute the sum of the integrals on the interface Si, We write = Jk, enu- 



k=l 



merating the terms in the same order of the list in ([7]). For each k = 1, . . . , 10 we denote 
as [Jfe] the sum of the /c-term coming from the integrations by parts in 0,i with the 
corresponding one from 0,2- 
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Proposition 5. If tp satisfies hypothesis (Hi), {H2) and (Tr) then there exist Aq and 
sq such that 
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Jsi = Y.^Jk{wi) + Jk{w2)) > 



(11) 



k=l 



for all V £ Z , yX > Aq, s > sq. 



Proof : 

It is not difficult to clieck tliat [Jk] = for k = 1,2 since ip and w satisfy tlie 
transmission conditions (Tr). Moreover, ip is constant on tlie interface and then we 

dip dw 

obtain Vip • Vw = -r— on Ti. Tlierefore, thanks to (.^2) we get 



-2sX 1 1 e 



ai 



dwi 



'^ + '^\dadt 



> s\5 





dwi 









dadt. 



dw 



By mean of the orthogonal decomposition Vw = -^r—v + VrW, where VrW is the 

ou 

projection of Vw on the tangent hyper-plane of dO,i, and from hypothesis (-^2) and 
(i^s) and the fact that Vt-wi = 'VrW2 we obtain 



sX 



ai 



dwi 



^J^ + '^Uladt 



sX 



dvi 

2 ( 2^^^! 2^V'2 \ 

IVrUil ( Oi";^ 1-02";^^ — I d(^di 



sX 



dv\ dv2 



V 5^ 1 dadt 



dui dv2 



_L \ / ; fliT7 |2 I dipi dip2 ^ , 
+ sX 0\\/rW\ aia2 I —— — 1 dadt 



> sX 

> 2sX5 



ai 



dwi 



V 5^ 1 dadt 



ai 



dvi 



dadt 



[Je] 



-s'X' 



ai 



dipi 



dvi 



^ + ^]dadt 

OUi OU2 



9^a\\'wi\^ dadt 



> s^X^6^ 

Since a G W'^'°°{Q) and £ C"'(ll), we also have 

mi < cs^x^ 



Iwil^e-'dadt + CX 



dwi 



>dadt 
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and 



10 



k=7 



wi\^9'-^ dadt. 



Thus, for s large enough, we get the desired result 

J2iJk]>{sS-c)(s^x' [[ \wi\^e^ + x 



dwi 



^ > 



2.2.4 Carrying all together. 

From ([To]) and Propositions [4] and [5] we obtain 



\w\\ — C \\w\\ — sXC 

s,X,ip s,X,ip,U 



dw 
ou 



<CRe{Pi{w),P2{w))L2 (12) 



C 



Adding -{\P^{w)\l^ + |P2Hli2) to both sides of ^ we obtain 



|(|PiHli. + |P2Hli.) + ||^|i;,,,, 



c\\wr -sXC 



dw 



<C\P{w)-qw\l,, 



what means that for all s > S2 and A > A2, since C > is a generic constant, 

\Pl{w)\l2 + \P2{w)\l2 + \\w\" 



<C\P{w)-qw\l,+C\\w\\l^^^^^+sXC 



dw 



(13) 



Now, if {ip^,ip'^) is an e-pair of transmission weight functions (see Definition [2D , we 
have an estimate like ()13p for each ip'' with U = B^{xk) where Xj G $7, j = 1,2 and 
e > 0. 

We sum up both estimates and we can show that the left hand side of each inequality 
can absorb the right hand side term || • \\sXip'' ,Be{xu) fro™ other inequality provided 
that £ is small and A is large enough. Indeed, by assumption we have that tp'^—ip^ > 6 > 
in Be{xi). 

Then, by taking A large enough we have 

gA(v.2-v.i) > 2C in Be{xi) 

i.e. 

Ce^ < in B,{xi) 

and we conclude that ""^ of the right and side is absorbed by the term Utt'^H^ j 

of the left hand side. It is clear that an analogous result is true by interchanging and 
■0^. Theorem [3] is proved. H 
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2.3 Particular case. 

In this part of the work, we construct exphcit weight functions adapted to particular 
discontinuous coefficients. 

We need the fohowing definition. 

Definition 6. We say that the open, bounded and convex set U C (N > 2) is 
strongly convex if dU is of class and all the principal curvatures are strictly 
positive functions on dU . 

Remark 1. Let us note that U C is strongly convex if and only if for all plane 
n C M.^ intersecting U, the curve H H dU has strictly positive curvature at each point. 
In particular, a strongly convex set is geometrically strictly convex. 

We assume that C $7 is a strongly convex domain with boundary Fi of class C^, 
and we set = \ Oi. Thus, we have = F U Fi, where this is a disjoint union. 
We deal with the case where a is locally constant 

with Oj > for j = 1, 2. 

In order to construct a convenient weight function, take xq G and for each x G 
Q, \ {xq} define i{xo,x) = {xq + X{x — xq) : A > 0}. Since Qi is convex there is 
exactly one point y{x) such that y{x) G Fi n i{xo,x). Thus, we can define the function 
p;Q\ {xo} M+ by: 

p{x) = \xo -y{x)\. (15) 

Let e > be such that Br C i^i (and small enough in a sense we will precise later) and 
let < ei < £2 < £• Then we consider a cut-off function rj G C°°(M^) such that 

< ?] < 1, 1] = in Be-^{xo), 1] = 1 in ^}\Be^{xo). 

For each j G {1,2} we take k such that {j,k} = {1,2} and we define the following 
functions in the whole domain 

Qfc |_ _ |2 

Xx)^ 

where Mi and M2 are positive numbers such that 

ai - 02 = Ml - M2. (16) 
Then, the weight function we will use in this work is 

{Tpiix) X G ill 
(17) 
Tp2[X) X G ^l2- 

Throughout the paper, we will use the notations a(x) = 02101(2;) + 011^2(2^) and 
M = Ml on Oi and M = M2 on Q2, so that we can write 

ijj(x) = r](x)a(x) ---5 h M. 

/3(xj^ 

As we can see in the following result, the main property of the weight function is a 
consequence of the strong convexity of the interior domain Oi. 



ipj{x)=r]{x) ^ |x-xoP + Mj X G il. 
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Lemma 7. If C is strongly convex and if the function — >R+ is defined 

by fJ-{x) = — — — then Z)^/i^(x) is positive definite for all x G M^\{a:;o}, uniformly in 

p[x) 

bounded subsets o/]R^\{xo}. 
Proof : 

We shall deduce this Lemma from well-known properties of compact convex subsets 
of (called convex bodies). However, for the sake of completeness of this paper, we 
include in the Appendix a self-contained proof of this result. 

Assuming without lost of generality that = 0, it is not difficult to see that /z is 
the gauge function of the convex set Qi (in other words, /i is a seminorm whose unit 
ball is Oi, see [26], p. 43, and section 2.3 of [TO]). 

The proof that /i is a convex function of class (hence D^/i > 0) can be read 
in [To] (Theorem 2.1). Moreover, it is proved that for each x G M^\{0} the only null 
eigenvalue of D'^fj,{x) corresponds to the direction x (which is the radial direction). 
The others eigenvalues, as functions of x, are bounded below by a positive constant, 
uniformly in x 7^ given in a bounded subset of M^. 

Thus, there exists 6 > such that for all x £ Q we have 

D'^fi{x){v,v) >5\v\^ yv£x^ = {y£W : yx = 0}. (18) 

1 X 1 

On the other hand, we have V/x(x) = - — - -|- |x|V— (x). Since p is constant in the 

p[x) \x\ p 

radial direction, we get x ■ V — (x) = 0. Hence we deduce that 

P 

X ■ Vpix) = = /u(x) 7^ 0. (19) 
p(x) 

Take x,f G M^\{0}. Then v = vi-^^ + V2y^ where y is an unitary element of x-*-. 

In view of the fact that D'^p'^{v,v) = 2pD^p{v,v) + 2\v ■ Vp\^, from ([18|), and 
L>^/i(x)(x, x) = we get 



D^p\x){v,v) > 2fi{x)D^p{x){v2y,V2y) + 2 

,2 

12 



X 

vi— ■ V/i(x) 

X 



2 



> 2p(x)6vi + 2-%,fi'^(x) 

xr 



2 

2fi{x)6vl + -^^vl 



and we conclude that D'^p?{x) is positive definite. 

■ 

Assuming the additional hypothesis about the sign of the jump on the interface, we 
can prove that the functions we have defined work as a weight function: 
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Proposition 8. Let Qi be an open and bounded set in M.'^ with smooth boundary, and 
fli, a2 real numbers such that: 

1. Vti is strongly convex. 

2. < 02 < ai- 

Then, for each pair of points xi,X2 G i^i, there exists e > such that the above con- 
struction gives up an e-pair of transmission weight functions {tfj^ in the sense of 
Definition [H 

Proof : 

For xo G f^i let ip be the function constructed as above and defined by (jl7p . If x € Ti 
we have p{x) = \x — xq\ and ^pj{x) = + Mj. From (fTUj) we get ipi = 'ip2 = c on Fi. 
Moreover, if x G Fi we have 

(I ^ — "^0 1 ^ \ 
—-^-^1 = a2V^p2{x). 

Hence (Tr) and (Hi) are satisfied (recall that z^i = —1^2 on. Ti). 

On the other hand, since Fi is a level set of ipi, then tpi{x) < 02 + Mi < ipi{y) for 
any x G Qi and y G O25 and we have — — > on Si and 



dvi dv2 dui \ a2 
what gives {H2)- 

For X G J7o\^e(^o)i denoting c{x) = ^^^^ , we get VV' = 2c(2;)(x— 3;o)+|x— 2;opVc(x). 
By construction c{x) is constant in the direction of x — xq, hence 

(x — xq) • Vc(x) = 

and then 

jV^Ap = 4c2(x)|x-xoP + |x-xohVc(x)p 
> 4c^(x)|x — xol^. 

Thus we have 

2 

2 



in r^o \ -Be(xo) and satisfies (H^) in that set. 
Property (i?4) is deduced from Lemma [71 

One can notice that xq can be arbitrarily chosen in since it is convex. Therefore, 
we can take two different points xi, X2 in Qi and we can construct the respective 
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weight functions il^^ and 'i/'^- For each k = 1,2, if)^ is a transmission weight function in 
\ -Be(xfc) and it remains to be shown that {H^) is fulfilled in order to finish the proof 
of Proposition [HI 

Let be d = — X2I such that e < d. On the one hand, for all x G Bi;{xi) we have: 

i^H^) < -^e^ + M < -^e^ + M, 
Pi "I 

where ai = d{xi,Ti) > 0. On the other hand, if we denote D2 = maxd{y, X2), we get, 
for all X G B^{xi), 

^'^(x) > ^d^ + M> -^d^ + M. 
P2 ^2 

Consequently, we have 

>af^-4l VxG5,(xi). (20) 

It is clear that an analogous result is true by interchanging xi and X2 (now with 02 and 
-Di). Thus, taking e < min we get (/fs) and Proposition [8] is proved. ■ 

^From Proposition [8] and Theorem [3l we obtain the following result: 

Theorem 9. Let the coefficient a be constant in the open set ilj and equal to aj for 
each j = 1,2. Suppose that 02 < ai and that Qi is an open, bounded and strongly 
convex set with smooth boundary. Then we have a Carleman estimate like Q for the 
Schrddinger equation ^ in the domain VL. 

3 Stability of the inverse problem 

As described in the introduction, will only give a local answer about the determination 
of the potential p. We will first work on a linearized version of the problem and consider 
the following Schrodinger equation : 

iu' + div (a(x)Vu) + q{x)u = f{x)R{x, t), 17 x (0, T) 
M(x,t) = 0, dQ.x{Q,T) (21) 

m(x,0) = 0, 

Here we set y = y{p) the weak solution to ([I]) and u = u{f) the one to ()21|) . If we 
formally linearize equation ([1]) around a non stationary solution, we obtain equation 
(pTI) . In fact, we notice here that if we set / = p — g, u = y{q) — y{p) and R = y{p) 
on X (0,T), we obtain (pT]) after substraction of ([H) with potential p from ([1]) with 
potential q and linearization. 

Linear inverse problem : Is it possible to determine f\Q from the knowledge of 
the normal derivative f^|gQx(o T) '^'^^'^^ P given and u is the solution to ([2T]) ? 

The following theorem proves that this inverse problem is well posed. 
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Theorem 10. Let q G L°°(J7) and u be a solution of equation l\21\i . We assume that 

R e W^'^{0,T,L'^{n)), 

i?(0) is real valued and \R{x,0)\ > > 0, a.e. in Q. 
There exists a constant C = C{Vt,T, \\q\\L°°[n)^R) > such that if 

du 



— (^H\0,T;L\ro)), 



then, 



Proof 



L2(r2) 



< C 



0.2 



du 
du 



H^{0,T;L^{dn)) 



(22) 



du . 



As we need to estimate — in H^(0, T; L'^CTq)) norm, we work on the equation satisfied 

01/ 

hy V = u' : 



iv' + div {a{x)Vv) + q{x)v = f{x)R'{x, t), n x (0, T) 
v{x,t) = 0, dnx{0,T) 
vlx,0) = -if{x)R{x,0), n 



(23) 



The Carleman inequahty we just obtained is the key of the proof. We extend the func- 
tion v on fl X {—T,T) by the formula v{x,t) = —v{x, —t) for every {x,t) £ Q x (— T, 0). 

(JV 

Since i?(0) and / are real valued, v £ C{[-T,T]; H^{n)) and — G L^{{-T,T) x T). 

We also extend R on Q x {—T,T) by the formula R{x,t) = R{x,—t) for every 
(x, t) € Q X (— T, 0) and if we denote the extention of R' by the same notation, then 
R' G L'^{-T,T;W^''^{n)). Thus, v satisfies the same equation 1^, set in {-T,T). 

As defined in Theorem [31 for A; = 1, 2, we set w = e~'^^ v and 

pfw^ = idtw^ + div(aVu;^) + s^alV/pw;^ 

Therefore, we define the following: 

2 



V/m / / Pi 



w dxdt. 
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On the one hand 

2 .0 



V'/m / [ pf w'' w'^ dxdt 



k=i 

2 



k=i -^-T 

2 

E 

k=l 



idtw^ + d\v{aVw^) + s^a\V (p^\^w^] dxdt 



2 .0 



Re l^dtw^w^ j — Im ( a 
dt dxdt 



2 .0 



fc=i 

2 



dxdt 



k=i •'^ 



(24) 



On the other hand, Cauchy-Schwarz inequahty and Carleman estimate from Theo- 
rem [9] give : 



/ < iz^f I ipfw'i^dxdty (^1' I iw^'i^dxdt 



< 



Eh' 



* 'w" ) 



k=l 



< Cs-l^[\\P^\w'' 



w 



LHQ) 



k=l 
2 



+ S 



dv 



dadt 



k=l 



< // \fR'\^e~^'^ dxdt + s 



dv 



a2 



dv 



dadt . 



Then, ^^{x,t) = ^T-t){T+t) '^^ ^^^^ e-2^^'(^'*) < g-^^^'C^.o) for all x G and 

t G {—T, T) and it is easy to see that 9e~'^^^ is bounded on E^ and that using the 
definition of the extensions of v and R' , we easily get 



02 



5u 



do-dt (25) 



^From R G VF^'''(0, T, L°°(0)) and |i?(x,0)| > ro > almost everywhere in Q, we 
deduce that 



3 <7o G ^'(0,T), |i?'(x,t)| < 5o(t)|i?(2;,0)|, Vx G f], t G (0,T). 
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Hence, from (p^ and (p5]) we have : 

2 



5^/ |/ni?(0)pe-2^^'=(°)dx 
13 -'n 



fc=i 

k=l \ £1 S 

< ^E' ^ 



dv 



dadt 



' ^ Tf \fW\R{Ore-'^^'''^'Uxdt 
Jo Jn 



k=l 





dv 







dadt. 



k=l 

But go £ -L2(0,r) implies /" \goit)\^ dt < K < +oo and 



so we can write 



s2 J ^^^Jn ^^,JJt 





dv 







dadt 



that becomes easily, if s is large enough {s > (C-fC) ^ ) and C remains a generic positive 
constant 

|/|2|i?(0)|2 (e-2V(0) + e-2V{0)\ < Cs-"2 [[ dadt. 
n ^ ' J Jt. ov 

Moreover, since |i?(x,0)| > ro > and e~ ' > e ^r^" > o almost everywhere in 

n, we obtain 



[ \f{x)\''dx<C 
Jn 

and therefore, Theorem [10] has been proved. 



^2 



dv 
dv 



dadt, 



Remark : if we replace the assumption "i?(0) is real valued " by the following 
"i?(0) takes its values in iM" , then the appropriate extensions for (x, t) in 0, x (— T, 0) 
are v{x, t) = v{x, —t) and R{x, t) = —R{x, —t). 

We will end this paper by the proof of Theorem [T] which is a direct consequence of 
Theorem [lOl Indeed, if we set u = y{q) — y{p), f = p — q and R = y{p), then u is the 
solution of 

u' + div(aVn) + (p - f)u = f{x)R{x, t) (0, T) x 
u = (0,r)xS (26) 

u(0) = n 

where q = p — f £ V(, with U bounded in L°°(J7) from the hypothesis of Theorem [TJ 
The key point is that in the proof of Theorem \T0\ all the constants C > depend on 
the L°°-norm of the potential. Thus, with q € U, we are actually, with equation (j26p in 
a situation similar to the linear inverse problem related to equation ()2ip and we then 
obtain the desired result. H 
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Appendix: Direct proof of Lemma [7] 

Without lost of generality, we can take xq = 0. Now, take x,v ^ \ {0} and define 
g{t) = fj?{x + tv) for t G M. Then g depends only on the restriction of /x^ to the plane 
n = {{x, v}) C spanned by the vectors x and v. Moreover, by definition of p, it is 
not difficult to see that pin = Po, where we have denoted by po the function defined in 
the plane 11 as in p^ . but where the closed curve is given by Ti = 11 n d^i, wich by 
hypothesis is strongly convex (see Remark [T|). 

It is not difficult to see that -r^(0) = D^{a^){x){v,v) and then this expression 

depends only on the curve Ti C H. We conclude that it suffices to consider the two- 
dimensional case. 

Assuming = 2, Fi can be parameterized in polar coordinates by 

7(0) = {p{e) cos e, p{e) sin e) ^ g [o, 2^) . 

The expression for the Hessian matrix of second derivatives in polar coordinates is 

D\p^) = QeH{p^)Qj 
where Qg is the rotation matrix by angle 6, and 



dr^ r \ drdd r 86 

\ r \drde r 89 J 7^~de^ r dr 



Now, since xq = 0, we have fi^(0, r) = + M. 

p{ey 

One can notice that /i^ is well defined and smooth in JIq \ ^^{xq) (which means {r > 
e} \ Ti). All the computations that follows are valid in this set. We already said above 
that p is constant with respect to r and only depends on 6 such that §f = 0. Hence, we 
get 



where we have denoted Pe = ^ and so on. 

We will use the following well known lemma (see [12]) concerning curves in the plane: 

Lemma 11. Let 7 6e a curve in the plane parameterized in polar coordinates by its 
angle: ^{6) = {r (9) cos 9, r (9) sin 9). Then, the curvature of ^ is given by the formula 

_ r^ + 2rj - rree 
(r2+r2)3/2 • 

Since the polar parametrization of Ti is given by r{9) = p{9) and VLi is strongly 
convex, we obtain 

(^) = (p2+ ^2)3/2 > ^ 
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We will now prove that H{fi'^) is uniformly positive definite in Qq \ B^{xq), which 
will imply that the hypothesis (-^4) is fulfilled in this set, since a is piecewise constant. 



2 

The eigenvalues of the matrix i^H^fi^) satisfy the equation 



r = — lb \J (fi — 4m 
where d = — (3pg — ppeg + 2/?^) and 

m = - {2pl - ppee + p") = t^'^t, (^) > 0. 

Then r2 = -(d + \/ d^ - 4m ] < d, and n = — > — for ah e [0, 27r). Since fii 
2 V / r2 d 

is bounded we get the desired result. H 
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